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1. Introduction 

In this paper we develop lower and upper bounds for 
arrowhead matrices. A matrix Q G IR mxm is called an 
arrowhead matrix if it has a form as follows: 

D c~ 


where D G a diagonal matrix, c is a vector 

in IR.” 1 ^ 1 , and b is a real number. Here the superscript 
“t” signifies the transpose. The arrowhead matrix Q is 
obtained by bordering the diagonal matrix D by the vector 
c and the real number b. Hence, sometimes the matrix 
Q in (1) is also called a symmetric bordered diagonal 
matrix. In physics, arrowhead matrices have been used to 
describe radiationless transitions in isolated molecules [1] 
and oscillators vibrationally coupled with a Fermi liquid [2]. 
Numerically efficient algorithms for computing eigenvalues 
and eigenvectors of arrowhead matrices were discussed in 
[3]. The properties of eigenvectors of arrowhead matrices 
were studied in [4], and as an application of their results, an 
alternative proof of Cauchy’s interlacing theorem was given 
there. The existence of arrowhead matrices was investigated 


recently in [5-8] such that the constructed arrowhead 
matrix has the pregiven eigenvalues and other additional 
requirements. 

Our motivation to study lower and upper bounds of 
arrowhead matrices is from Kung and Suter’s recent work on 
the hub matrix theory [9] and its applications to multiple- 
input and multiple output (MIMO) wireless communication 
systems. A matrix, say A, is a hub matrix with m columns if its 
first m— 1 columns (called nonhub columns) are orthogonal 
to each other with respect to the Euclidean inner product 
and its last column (called hub column) has a Euclidean 
norm greater than any other columns. Subsequently, it was 
shown that the Gram matrix of A, that is, Q = A*A, is 
an arrowhead matrix and its eigenvalues could be bounded 
by the norms of the columns of A. As pointed out in 
[9-11], the eigenstructure of Q determines the properties 
of wireless communication systems. This motivates us to 
reexamines these bounds of the eigenvalues of Q and makes 
them sharper. In [9], the hub matrix theory is also applied 
to the MIMO beamforming problem by comparing k of m 
transmitting antennas with the largest signal-to-noise ratio, 
including the special case where k = 1 which corresponds 
to a transmitting hub. The relative performance of resulting 
system can be expressed as the ratio of the largest eigenvalue 
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of the truncated Q matrix to the largest eigenvalue of the 
Q matrix. Again, it was previously shown that these ratios 
could be bounded by the ratios of norms of columns of the 
associated hub matrix. Sharper bounds will be presented in 
Section 4. 

The well-known result on the eigenvalues of arrowhead 
matrices is the Cauchy interlacing theorem for Hermitian 
matrices [12]. We assume that the diagonal elements dj, 
j = 1,2,..., m — 1, of the diagonal matrix D in (1) satisfy 
the relation d\ < d 2 < ■ ■ ■ < d m -i- Let Ai, A 2 ,..., A m be the 
eigenvalues of Q arranged in increasing order. The Cauchy 
interlacing theorem says that 

^1 — d\ < A 2 ^ — ■ ■ ■ — d m —2 — A m _i < d m -\ < A m . (2) 


Lemma 1. Let Q E K. mxm be an arrowhead matrix of the form 
(1), where D = diag(di,d 2 »--->^m-i) £ ||( m -i)x(m-i) j b E R, 
andc = (ci,C2,...jC m -i) E [R m_1 . Then 

m— 1 m -1 2 m— 1 

det(Xl-Q) = (X-b)Y\(X-d k )- X \cj[Y\( X - d kl 

k =1 ;' = 1 k=l 

kfj 

( 6 ) 

The proof of this result can be found in [5, 13] and 
therefore is omitted here. 

When the diagonal matrix D in (1) is a zero matrix, the 
following result is followed from Lemma 1. 


When the vector c and the real number b in (1) are taken 
into consideration, a lower bound of Aj and an upper bound 
of A m were developed by using the well-known Gershgorin 
theorem (see, e.g., [3, 12]), that is, 


m— 1 


A m < max|iii + |ci d m -i + \c m -i\,b+ X kilj, (3) 

Ai > mint d\ - |ci |,..., d m -i - |c m _i|,fc- X Id \- ( 4 ) 


i= 1 


Accurate bounds of eigenvalues of arrowhead matrices 
are of great interest in applications as mentioned before. 
The main results of this paper are presented in Theorems 
11 and 12 for the upper and lower bounds of the arrowhead 
matrices. It is also shown in Corollary 13 that the resulting 
bounds are tighter than in (2), (3), and (4). 

The rest of the paper is outlined as follows. In Section 2, 
we will introduce notation and present several useful results 
on the eigenvalues of arrowhead matrices. We give our 
main results in Section 3. In Section 4, we revisit the lower 
and upper bounds of the ratio of eigenvalues of arrowhead 
matrices associated with hub matrices and wireless com¬ 
munication systems [9], and subsequently, we make those 
bounds shaper by using the results in Section 3. In Section 5, 
we compute the bounds of arrowhead matrices using the 
developed theorems via three examples. Conclusions are 
given in Section 6. 


Corollary 2. Let Q E M. mxm be an arrowhead matrix having 
the following form: 


Q 



(7) 


where c is a vector in IR m 1 and b is a real number. Then the 
eigenvalues of Q are 


b - ^b 2 + 4||c|| 2 b + -\jb 2 + 4||c|| 2 

2 > — 2 > 

A;(Q) = 0, for i = 2,..., m - 1. 

( 8 ) 


Proof. By using Lemma 1, we have 

det(A7 - Q) = A m ~ 2 (A 2 -bX- ||c|| 2 ). (9) 

Clearly, A = 0 is a zero of det(A/ - Q) with multiplicity m — 2. 
The zeros of the quadratic polynomial A 2 - bX - ||c|| 2 are 
(b - fb 2 + 4||c|| 2 )/2 and ( b + fb 2 + 4||c|| 2 )/2, respectively. 
This completes the proof. D 


In what follows, a matrix Q having a form in (7) is called 
a special arrowhead matrix. The following corollary (also, see 
[3]) is a direct result from Lemma 1. 


2. Notation and Basic Results 

The identity matrix is denoted by I. The notation 
diag(fli, fl 2 > • • • > (in) represents a diagonal matrix whose diag¬ 
onal elements are a\ , < 22 ,..., a n . The determinant of a matrix 
A is denoted by det(A). The eigenvalues of a symmetric 
matrix A E R nXn are always ordered such that 

Ai (A) < A 2 (A) < ■ ■ ■ < A „{A). (5) 

For a vector a E K”, its Euclidean norm is defined to be 
Ha|| := I a,-1 2 . 

The first result is about the determinant of an arrowhead 
matrix and is stated as follows. 


Corollary 3. Let Q be an m X m arrowhead matrix given by 
(1), where D = diag(di,d 2 ,---.d m -i) E [rC 1 ”—i)x(m— 1 )^ b E K, 
and c = (ci,c 2 ,..., c m - 1 ) E IR m_I . Let us denote the repetition 
of the number dj in the sequence {d,-}” j by kj. If kj > 2, then 
dj is the eigenvalue of Q with multiplicity kj — 1. 

Proof. When the integer kj > 2, the result follows from 
Lemma 1 since (A — dj) ki l is a factor of the polynomial 
det(A7 - Q). Id 

Corollary 4. Let Q be an m x m arrowhead matrix given 
by { 1), where D = diag(di, d 2 ,... ,d m -i) £ |j^(m-i)x(m-i) ) 
b E M, and c = (ci, C 2 ,..., c m - 1 ) e r^ 1 . Suppose that the 
last k > 2 diagonal elements d m _^, d m -/t+ 1 , ■ ■ ■ > d m -1 of D are 
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identical and distinct from the first m-k - l diagonal elements 
d\,di,..., d„,-it_i ofD. Define a new matrix 


[di 


Cl 


Q:= 


Cl 


d-m-k-l Cm—k-1 

d-tn—k C m -k 
Cm—k —1 Cm-k b 


( 10 ) 


with Cj = Cj for j = 1,2,... ,m — k — 1 and c m - k = 
■\jllf=m-k \ c j\ 2 - Then the eigenvalues ofQ are that ofQ together 
with d m -k with multiplicity k - 1. 


entries of the vector c in (1) are nonzero. The reason for this 
assumption is the following. Suppose that C;, the jth entry of 
c, is nonzero, it can be easily seen from Lemma 1 that A - dj 
is a factor of det(A/ - Q); that is, dj is one of eigenvalues of 
Q. The remaining eigenvalues of Q are the same as those of 
a matrix which is obtained by simply deleting the jth row 
and column of Q. In summary, for any arrowhead matrix, 
we can find eigenvalues corresponding to repeated values in 
D or associated with zero elements in c by inspection. 

In this paper, we call a matrix Q in (1) irreducible if 
the diagonal elements di, d 2 ,..., d m -i of Q are distinct and 
all elements of c are nonzero. By using Corollary 4 and the 
above discussion, this arrowhead matrix can be reduced to 
an irreducible one. 


Proof. Since numbers d m -ib d m -/t+i,... , d„,-i are identical 
and distinct from numbers di, d 2 ,... ,d m - k - 1 > we have 


m-1 


n (^ - di) 

i=i 


/ \ 

m-k 

n(A- di) 
1=1 

Vti ) 


(A -d m . k ) k ~\ 


j <m-k - 1, 


Remark 5. In [4, 9], Hermitian arrowhead matrices are 
considered; that is, it allows that c in the matrix Q of the form 
(1) is a vector in We can directly construct many (real 
symmetric) arrowhead matrices denoted by Q from Q. The 
diagonal elements of these symmetric arrowhead matrices 
are the exactly same as those of Q. The vector c in Q could 
be chosen as 


m-1 2 m 

I |Cj I UiX-di) 

j=m—k i=l 

i*j 

( m-1 2 \ / m-k—1 \ 

x | Cj | n (A-d,o \(x-d m - k ) k -\ 

j=m—k J \ i= 1 / 

(ID 


c= (±|Ci|,±|C 2 l,...,±|c m -i|). (13) 

In such a way, there are 2 m ~' such symmetric arrowhead 
matrices. Because det(A7 - Q) = det(Af - Q) by Lemma 1, 
every such symmetric arrowhead matrix Q has the identical 
eigenvalues with Q. This is the reason why we just consider 
the eigenvalues of real arrowhead matrices in this paper. 


By (6) in Lemma 1, we have 
det(A7 - Q) 


The following well-known result by Weyl on eigenvalues 
of a sum of two symmetric matrices is used in the proof of 
our main theorem. 


m-1 m-1 2 m_ ^ 

= (\-b)Y\tt-di)- X \ c i | m-di) 

i= 1 ; = 1 i= 1 

it) 


Theorem 6 (Weyl). Let F and G be two m x m symmetric 
matrices. Let us assume that the eigenvalues of F, G, and F + G 
have been arranged in increasing order. Then 


( 

m-k 

(x-b)Yl(x-di)- 

i=i 

V 


m-k 

1 


;=i 


.m-k-l 

n 


\ 

(A -dd 


2=1 

/ 


(A -d m - k ) k ~ l 


= det(A7 - Q) ■ (A -d m -k) k ~\ 


( 12 ) 


Clearly, if A is an eigenvalue of Q, then A is either an 
eigenvalue of Q or d m - k - Conversely, d m - k is an eigenvalue 
of Q with multiplicity k - 1 and the eigenvalues of Q are that 
of Q. This completes the proof. EUl 

By using Corollaries 3 and 4, to study the eigenvalues of 
Q, we may assume that the diagonal elements d\, d 2 , ■ ■ ■, d m -i 
of Q are distinct when we study the eigenvalues of Q in 
(1). Since eigenvalues of square matrices are invariant under 
similarity transformations, we can without loss of generality 
arrange the diagonal elements to be ordered so that d\ < 
da < ■ ■ ■ < d m - 1 . Furthermore, we may assume that all 


Xj(F + G) < A,(F) + Aj_ i+m (G), 

for i> j. 

(14) 

A j(F + G) > Xj(F) + A ; -j + i(G), 

for i < j. 

(15) 

Proof. See [14, page 62] or [12, page 184] 


□ 


To apply Theorem 6 for estimating eigenvalues of an 
irreducible arrowhead matrix Q, we need to decompose Q 
into a sum of two symmetric matrices whose eigenvalues are 
relatively easy to be computed. Motivated by the structure 
of the arrowhead matrix and the eigenstructure of a special 
arrowhead matrix (see, Corollary 2), we write Q into a sum 
of a diagonal matrix and a special arrowhead matrix. 

To be more precisely, let Q e [R mxm be an irreducible 
arrowhead matrix as follows: 


Q 


D 


( 16 ) 
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4 

where d m e RL, D = diag(di, d 2 ,..., d m _i) with 0 < d\ < 
d 2 < ■ ■ ■ < d m _1 < d m , and c is a vector in [R m_1 . For a given 
p e [0,1], we write 

Q = £ + S, (17) 


Therefore, we can use Theorem 6 to give estimates of the 
eigenvalues of Q via those of E and S. To number the 
eigenvalues of E, we introduce the following definition. 

Definition 7. For a number p e [0,1], we define an operator 
T p that maps a sequence {d ; }J =l satisfying 0 < d\ < d 2 < 

■ ■ ■ < d m _1 < d m to a new sequence {dj}f =l := T p ({di}"‘ =1 ) 
according to the following rules: if pd m < d\, then di := pd m 
and dj+i := dj for j = 1 ,...,m — 1; if pd m > d m _i, then 
dj := dj for j = — 1 and d m := pd m ; otherwise, 

there exists an integer j 0 such that dj 0 < pd m < dj 0+ 1 , then 

dj . dj for j 1,..., j 0 , dj 0 +\ . pd m , and dj+i . dj for 

= j 0 + - 1. 

Theorem 8. Let Q e IR mxm be an irreducible arrow¬ 
head matrix having a form of (16), where D = 
diag(di,d 2 >-■ ■ >d m _i) with 0 < d\ < d 2 < ■ ■ ■ < d m -i < d m , 
and c is a vector in IR m_1 . For a given p E [0,1], define 
{di} n j l =l := r p ({d,} ; m =1 ). Then, one has 



min 

di 

+ t, d 2 , d m + sj, 

if j = 1. 



Aj(Q) < - 

min 

dj 

+1, dj+i j, 

if 2 < j < m 

- 1, 


d m + 

t, 


if j = m, 









(19) 


d\ + 

s, 


ifi = 

1, 


O 

IV 

max 

[dj 

-l, dj + sj, 

if 2< 

j * 

m - 1, 


max 

d\ 

+ t , d m _ i , d m I - s 

> fj = 

m, 



(1 - p)d m - ^(l-p) 2 d 2 m +4\\c\\ 2 

s = - - -, 

2 

( 21 ) 

= (l-p)d m +V(l-p) 2 ^+4||c|| 2 

2 

Proof. For a given number p e [0,1], we split the matrix Q 
into a sum of a diagonal matrix E and a special arrowhead 
matrix S according to (17), where E and S are defined by (18). 
Clearly, we know that 

A j(E) = dj (22) 
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for j = 1,2,..., m. By Corollary 2, we have 

Ai(S)=s, A m (S) = t, Aj(S) = 0, for j = 2,...,m - 1, 

(23) 

where s and t are given by (21). 

Upper Bounds. By (14) in Theorem 6, we have 
Xj(Q) < Xi(E) + X 

m+j- i(S ) (24) 

for all i > j. Clearly, for a given j, 

X ; (Q) < minjAiCE) +A m+j—i (S)}- (25) 

More precisely, since {d/}™ j is monotonically increasing, s < 
0, and t > 0, we have 

Ai(Q) < minjdi + t, d 2 ,..., d m + sj 

= min jdi + t, d 2 , d m + s|, 

A ; (Q) < min|dj + f,dj + i,...,d m } = minjdj + t,d j+ 1 | 

(26) 

for j = 2,..., m - 1, and 

A m (Q) ^ A m (E) + X m (S) = d m + t. (27) 

In conclusion, (19) holds. 

Lower Bounds. By (15) in Theorem 6, we have, for a given j, 
Aj(Q) > max|A,(E) + Aj_ ;+ i(S)|. (28) 

Hence, 

Ai(Q) > Ai(£) + Ai(S) = d\ +s, 

Aj(Q) > maxjdj + s,dj~ i,... ,di | = maxjdj + s,dj_i| 

(29) 

for j = 2,..., m — 1, and 

A m (Q) > maxjd m + s, d m -i,..., d 2 , d\ + t| 

= max|d m + s, d„,_i, d\ + f j. 

As we can see from Theorem 8, the lower and upper 
bounds of the eigenvalues for Q are functions of p e [0,1 ] for 
the given irreducible matrix Q. In other words, the bounds 
of eigenvalues vary with the number p. Particularly, when we 
choose p being the ending points, that is, p = 0 and p = 1, 
we can give an alternative proof of interlacing eigenvalues 
theorem for arrowhead matrices (see, e.g., [12, page 186]). 
This theorem is stated as follows. 
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Theorem 9 (Interlacing eigenvalues theorem). Let Q e R mXm 
be an irreducible arrowhead matrix having a form in (16), 
where D = diag(di, d 2 ,... ,d m -\) with 0 < d\ < da < ■ ■ ■ < 
d m - i < d m , and c is a vector in Let the eigenvalues of Q 
be denoted by withX\ < X 2 < ■ ■ ■ < X m . Then 

Ai < <4 < A 2 — d 2 ^ ^ d m —2 — A m _j < d m -\ < A m . 

(31) 


Prop/! In Theorem 8, the upper bounds of the eigenvalues of 
Q in (19) are determined by dj, j = 1,2,... ,m, and s and t 
in (21). They can be viewed as functions of p in [0,1], That 
is, the upper bounds of the eigenvalues of Q are functions of 
p in the interval [0,1]. Therefore, we are able to find optimal 
bounds of the eigenvalues of Q by choosing proper p. The 
upper bounds on XjiQ) for j = 1, 2 < j < m - 1, and j = m 
in (34) are discussed separately. 


Proof. By using (19) with p = 0 in Theorem 8, we have 
A j < dj for j = 1,2,..., m - 1. By using (20) with p = 1 
in Theorem 8, we obtain A j > dj- 1 for j = 2,3 ,...,m. 
Combining these two parts together yields our result. □ 

The proof of the above result shows that we could 
have improved lower and upper bounds for each eigenvalue 
of an irreducible arrowhead matrix by finding an optimal 
parameter p in [0,1]. Our main results will be given in the 
next section. 

3. Main Results 

Associated with the arrowhead matrix Q in Theorem 8, we 
define four functions f, i = 1,2,3,4, on the interval [0,1] as 
follow: 

flip) :=\((l-p)d m -^(l-p) 2 d] n +4\\c\\ 2 ), 

flip) := \((l-p)d m + ^l-p) 2 dl+m\ 2 ), (32) 

hip) := pd m +flip), 
hip) ■■= pd m + flip)- 

Obviously, 

s = fiip), t = f 2 ip), (33) 


Upper Bound ofXiiQ). From (19), we have 

Ai(S) < minjdi + f,d 2 ,(i m + sj, (35) 

where <4, s, and t are functions of p on the interval 
[0,1]. In this case, we consider p in the following four 
subintervals. [0, d\!dy^\, [d]/, di/d J, [ di/d ^^, d^ .\l d,,,], 
and [d m -i/d m , 1], respectively. For p e [0,di/d m ], we have 
d] + t = hip), di = di, and d m + s = flip). For p e 
[d]/d m ,di!d m ), we have d\ + t = d\ + flip), di = pd m , and 
d m + s = d m -1 + flip). For p e [d 2 /d m ,d m -\ld m ], we have 
di+t=di+ flip), d 2 = d 2 , and d m + s = d m -1 + flip). For 
p e [d m -i/d m , 1], we have d\ + t = di + flip), di = d 2 , and 
d m + s = hip). Hence 


min d 2 = di, 

pe[ 0,1] 



where s and t are given by (21). 

The following observation about monotonicity of func¬ 
tions f, i = 1,2,3,4, is simple, but quite useful as we will see 
in the proof of our main results. 

Lemma 10. The functions /] and f 2 both are decreasing while 
h and h ar £ increasing on the interval [0,1]. 

The proof of this lemma is omitted. 

Theorem 11. Let Q be an irreducible arrowhead matrix 
defined by (16) and satisfying all assumptions in Theorem 8. 
Then the eigenvalues of Q are bounded above by 


XjiQ) < 


minjdi.dm-! + /i(^ J ‘)}> 
dj, 


d m -1 + 



if j = 1, 

if 2 < j < m — 1, 


min(di + t 

pev ' 


fiiO), ifV = 

/ /L \ 

if V = 

dm— 1 


dl+f 2 it )’ 


0 

d\ d 2 
d, 


rfl + / 2 (^f)’ ifV “ 

di + fiil), if V = 


dy] 

'dm r 

dy\ 

'd m J’ 

dm- 1 "] 

> —j > 

Clm J 


(36) 


Since 0 > fi(di/d m ) > fiid 2 /d m ) > fi{d m -i/d m ), / 2 (0) > 
d m > d u and f 2 id 2 /d m ) > fiid m -i/d m ) > / 2 ( 1) > 0, we have 

Ai(Q) < mm^d u dm-i ( 37 ) 


if j = m. 


(34) 


5 


Upper Bound of XjiQ), for 2 < j < m- 1. From (19), we have 
Xj(Q) < minjdj + t,d j+ 1|. (38) 
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In this case, we consider p lying in the following four subin¬ 
tervals. [0, dj—\/djyf ], [dj—]/ dm, cij/dm |, [dj/dm , dj+i/dm], and 
[dj + }/d m , 1], respectively. For p e [0,dj-i/d m ], we have dj + 
t = dj-\ + f 2 {p) and dj+i = dj. For p e [dj-i/d m ,dj/d m ], we 
have dj + t = / 4 (p) and dj +1 = dj. For p e [dj/d m , dj+\/d m ], 
we have dj + t = dj + flip) and d, +1 = pd m . For p e 

[dj+i/d m , 1], we have dj + t = dj + flip) and dj+i = dj+\. 
Hence 


min dj+i = dj, 
pe[o,i] 


min(d i + t) = 


dj -i + /2 
dj-1 + /2 
dj + /2 


dj—i 
dm 
dj -1 


2j+l 


+ /2(1), 


if V = 


if V = 


ifV = 


ifV = 


0, 


dj-i 


dj- 1 dj 

d m d m 

dj dj+i 
dm d m 


dj -1 


,1 


(39) 


Therefore, 


jdj-i + / 2 ^ 


dj -1 


Aj(Q) < min-j dj- 1 + / 2 ( -y— ], d ; - [. 

Since d ; _i + / 2 (dj_r 
we get 

Aj(Q) < d ; -. 

Upper Bound of\ m {Q). From (19) we have 
Am ( Q) — ^m 3* t. 


(41) 


(42) 


For p e [0, d„,_i/d m ], we have d m + t = d m ~i + flip) while 
for p e [d m _i/d m , 1], we have d m + t = / 4 (p): 


min (d m + f) = 


peV 


Hence, 


dm-! +/ 2 r' 
V d m 

j , r l dm -1 


II 

> 

i i 

© 

3 i 

i_i 

II 

i i 

a. f»* 

3 7 

I— 1 

-T" l_l 

dm -1 \ 

dm / 

(' 


This completes the proof. 


□ 


Theorem 12. Let Q be an irreducible arrowhead matrix 
defined by (16) and satisfying all assumptions in Theorem 8. 
Then the eigenvalues of Q are bounded below by 


A AQ)> 


max-jdj-i, dj + fi{^~ 

dl+f 2 {t)’ 


if j = 1, 
if 2 < j < m — 
if j = m. 

(45) 


Proof. In Theorem 8, the lower bounds of the eigenvalues 
of Q in (20) are determined by dj, j = 1,2 ,...,m, and s 
and t in (21). As we did in Theorem 12, the lower bounds 
of the eigenvalues of Q are functions of p in the interval 
[0,1]. Therefore, we are able to find optimal bounds of the 
eigenvalues of Q by choosing proper p. The discussion is 
given for j = 1, 2 < j < m - 1, and j = m in (45), separately. 


Lower Bound o/Ai(Q). From (20), we have 
Ai (Q) > d\ + s. 


(46) 


In this case, we consider p lying in the following two 
subintervals: [0, d\/d m \ and [d\/d m , 1]. For p e [0,di/d m ], 
di +s = flip). For p e [dj/d m , 1], we have di -l-s = d\ + flip). 
Hence 


(47) 


(40) 

maxfdi + s) = - 

di + fi 

'H 

Him ' 

av = 

k j-l 

Um J 

— d m > dj , 

peV 1 > 

di + fi 

K Him ' 

av = 

ViA 

- HL m J 


It leads to 




Lower Bound ofXiiQ). From (20), we have 
A 2 (Q) > max|di,d 2 +s|. 


(48) 


(49) 


In this case, we consider p lying in the following three 
subintervals: [0,di/d m ], [di/d m ,d 2 /d m ], and [d 2 /d m ,l]. For 
p e [0,di/d m ], we have dj = pd m , d 2 + s = di + flip). For 
p e [di/d m ,d 2 /d m ], we have di = d\, d 2 + s = flip). For 
p e [di/d m , 1], we have d\ = d 2 and d 2 + s = d 2 + flip). 
Hence, 


max(d 2 + s) = 

peV ' ' 


maxdi = di, 

pev 


di + fdO), 
di 


Hi)- ■ 

Hi)’ ‘ 


d 2 + fi ( ~r~ 
di + fi ( -f- 


HV = 

if V = 

HV = 


0, dl\ 

d\ d-2 ~| 
dm d m J 

f.il. 

u m J 


( 50 ) 
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These lead to 

A 2 (Q) > maxjdi,d 2 + fi |- (51) 

Lower Bound ofXj(Q), 3 < j < m - 1. From (20), we have 
Aj(Q) > max|dj_i,dj + s|. (52) 

In this case, we consider p lying in the following three subin¬ 
tervals. [ 0 , dj— 2 /dm[, [dj— 2 / d m , dj — 1 / d m \y [dj ]/dm, dj/d m J, 
and [dj/d m , 1]. For p G [0, dj- 2 /d m ], we have dj_i = dj _ 2 
and dj 3- s = dj -1 + /i(p). For p G [dj- 2 /d m , dj-\/d m ], 
we have dj_i = pd m and dj + s = dj-\ + fi{p). For p G 
[dj_i/d m ,dj/d m ],wehavedj_i = dj -1 and dj+s = flip). For 

p G [dj/d m , 1], we have dj -1 = dj_i and dj 3- s = dj 3- /i(p). 
Hence 


maxd,_i = dj-i, 


max (dj 3- 5 ) = 


dj -1 3 - MO), 

if y = 

o, a -p 

dm 

j 



di—2 

di —1 

if y = 

j a 

d m 

j _ 

d m 

( dj\ 

\ u m / 


r , 

di 1 

if y = 

j j 

d m d m _ 

( dj \ 

+ /'.j ’ 

\ d m / 


r dj 1 


if y = 

dm 



(53) 


Since dj- 1 > dj_i + /i(0) > dj_i + /i(dj_ 2 /d„,), we have 


Aj(Q) > max]dj„i,dj + fi I '-j- 

V d n 


(54) 


Lower Bound of\ m (Q). From (20), we have 

A m (Q) > maxjdi + t,d m -\,d m 3-sj-. (55) 

In this case, we consider p lying in the following three subin¬ 
tervals. [ 0 ,di/d m ], [d| /d m , d m 2 /d m ], [d m — 2 /d m ,d m —i/d m \, 
and [d„,_i/d m ,l]. For p G [0,di/d m ], we have d\ + t = 
f^ip), d m -1 d m — 2 > d m 3" s d m -\ 3- f\ip). For p G 

[di/ dm, d m — 2 /dm ] , we have di 3- 1 di -t (p), dm —1 d m — 2 , 
d m 3- s — dm —1 3~ ftip)- For p G [dm— 2 /dm, dm— i/dm] , we have 
di 3- f = di 3- / 2 (p), d m -i = pd m , d m + s = d m -1 3- /i(p). For 


p £ [dm— 1 / d m , 1], we have di 3 -1 — d\ 3- flip ), d m -i — d m -i, 
d m 3 - s = flip). Hence 


max dm —1 — dm— 1 , 

pHO.l] 


max(d m 3- s) 

p£V v > 


max(di 3- f) 
pev v ’ 


d m -i + MO), 


d m -i+ fi ( -j- )> 


7 . r ( dm-2 

“m-1+ /l I “j- 

\ w m 

dm + /l(l), 

dl+f 2 it )’ 

di+/ 2 ( dl 
di +■ /a ^ 


d 

d m -2 


d\ 3- 


r (dm—l 

\ a m 


ifV = 

if V = 

, ifV = 

ifV = 

if y = 

if y = 

if y = 
if y = 


dm J 

d\ d m —2 
d m d m 

dm-2 


\ 


dm -1 1 

u m J 


— ll 

7 > 1 > 

w w J 

o,Al 

dm J 

^1 d m —2 


dm dm 

dm-2 d m 
dm dy, 


f]' 


-4 

1 J 


(56) 


Since 0 > /i(0) > /i(di/d m ) > /i(d m - 2 /d m ) and/ 2 (di/d m ) > 
fi(d m Jd m ) > / 2 (d m -i/d m ), we have 

A m (Q) > max|d m -i,d m 3 -/i(l),di 3 -/ 2 (^-) j. (57) 


Since di 3- / 2 (di/d m ) = Md\/d m ) > MO) > d m , we get 

Am(Q)>d 1+ / 2 (^). (58) 

This completes the proof. □ 

Corollary 13. Let Q be an irreducible arrowhead matrix 
defined by (16) and satisfying all assumption in Theorem 8. 
Then upper and lower bounds of the eigenvalues of Q obtained 
by Theorems 11 and 12 are tighter than those given by (2), (3), 
and (4). 


Proof. Since 


min|di,d m -i 3 -/i^^-d^| < di, (59) 


then the upper bound for the eigenvalue Ai(Q) given by (34) 
in Theorem 11 is tighter than that by ( 2 ). The upper bounds 
for the eigenvalues A jiQ), j = 2, ..., m - 1 , provided by (34) 
in Theorem 11 are the same as those by (2). 
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Note that 0 < d\ < ■ ■ ■ < d m -\ < d m ; the right-hand side 
of (3) with b = d m is 


I d\ + |Ci |, d m ~ i + \c m ~i I, b + X kil 


m -1 

= d m + I |. 
2=1 


(60) 


Since ||c|| < YULi kil, d m + ||c|l = Ml), and d m - X + 
fi{dm-i/d m ) = Mdm-i/dm), we have 


m -1 

dm I Cj | 
2=1 





(61) 


and then the upper bound of A„,(Q) from (34) in Theorem 11 
is tighter than that from (3). 

Now we turn to the lower bounds of Aj(Q). Since 




(62) 


4. Hub Matrices 

Using the improved upper and lower bounds for the 
arrowhead matrix, we will now examine their applications to 
hub matrices and MIMO wireless communication systems. 
The concept of the hub matrix was introduced in the context 
of wireless communications by Kung and Suter in [9] and it 
is reexamined here. 

Definition 17. A matrix A £ K" xm is called a hub matrix, if its 
first m— 1 columns (called nonhub columns) are orthogonal 
to each other with respect to the Euclidean inner product 
and its last column (called hub column) has its Euclidean 
norm greater than or equal to that of any other columns. We 
assume that all columns of A are nonzeros vectors. 

We denote the columns of a hub matrix A by 
«i, fl 2 , ■ ■ ■, flm- Vectors a\, ai,. ■. , a m -i are orthogonal to each 
other. We further assume that 0 < ||flil| < ||a 2 II < ■ ■ ■ < 
||a m ||. In such case, we call A an ordered hub matrix. Our 
interest is to study the eigenvalues of Q = A‘A, the Gram 
matrix A. In the context of wireless communication systems, 
Q is also called the system matrix. The matrix Q has a form 
as follows: 


for j = 2 ,..., m — 1 and 


d\ + 



^ dm dyi 


(63) 


llfltll 2 


\\a 2 \\ 2 


Q = 


{d\ , CL m ) 
(#2, tt-m ) 


we know that the lower bounds for the eigenvalues Aj(Q), 
j = 2,provided by (45) in Theorem 12 are tighter 
than those by (2). Ip 


II dm— 1II 

{fl-mi d\) {dm> CI 2 ) ’ " * {d m , Am— 1) 


1, ) 

II Q-m II 

( 66 ) 


Remark 14. When c in (16) is a zero vector, by using 
Theorems 11 and 12, we have dj < A(Q) < dj, that is, 
A(Q) = dj. In this sense, the lower and upper bounds given 
in Theorems 11 and 12 are sharp. 

Remark 15. When Q in Theorems 11 and 12 has size of 
2 x2, the upper and lower bounds of its each eigenvalue are 
identical. Actually, from Theorems 11 and 12 we have 


Clearly, Q is an arrowhead matrix associated with A. 

An important way to characterize properties of Q is in 
terms of ratios of its successive eigenvalues. To this end, the 
ratios are called eigengap of Q which are defined [9] to be 


EGi(Q) 


Am-(i-l) (Q) 
A m -i(Q) 


(67) 


d _1 + /1 - MQ) ^ minjdi,di + /1 ( . ) j> 

(64) 

dl+/2 (f) ^A2 ( Q)<di+/2(|). 

Clearly, we have 

Ai(Q) = di+/i(|), A 2 (Q) = d 1 +/ 2 (|). (65) 

This can be verified by calculating the eigenvalues Q directly. 


for i = 1,2 ,...,m — 1. Following the definition in [9], we 
define the ith hub-gap of A as follows: 

HG,(A) = (68) 

I a m -i II 

for i = 1,2 ,..., m — 1. 

The hub-gaps of A will allow us to predict the eigen- 
structure of Q. It was shown in [9] that the lower and upper 
bounds of EGi (Q) [9] are given by the following: 


Remark 16. For the lower bound of the smallest eigenvalue HGi(A) < EGi(Q) < (HGi(A) + 1)HG2(A). (69) 

of an arrowhead matrix, no conclusion can be made for the 
tightness of the bounds by using (4) and (45) in Theorem 12. 

An example will be given later (see Example 22 in Section 5). 
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These bounds only involve nonhub columns having the two 
largest Euclidean norms and the hub column of A. Using 
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the results in Theorems 11 and 12, we obtain the following 
bounds: 


/ 4 (||ai|| 2 /||fl„,|| 2 ) 


I ttm-1 I 


< EGi(Q) < 


M 

(II ^m- 1 

llVlkJI 2 ) 

maxj 

| II &m— 

2II 2 > / 3 I 

[\\a m M\ 2 /\\a m \\ 2 )} 


(70) 


Obviously, these bounds are not only related to two nonhub 
columns with the largest Euclidean norms and the hub 
column of A but also related to the nonhub column having 
the smallest Euclidean norm and interrelationship between 
all nonhub columns and the hub column of A. As we 
expected, the lower and upper bounds of EGi(Q) in (70) 
should be tighter than those in (69). To prove this statement, 
we give the following lemma first. 

Lemma 18. Let cq, a 2 ,..., a m be the columns of a hub matrix 
A with 0 < ||fli|| < ||a 2 ll < ■ ■ ■ < ||fl m -ill ^ ll«mll- Then 


Mp) > llflmll 2 for p e (0,1], 


< llflmll 2 + lla„,-il |2 


(71) 

(72) 


Proof. From Lemma 10, we know, for p E (0,1], 

,,v lkml| 2 +Vll«J| 4 +4||c|| 2 .. ||2 (73) 

Mp) > Mo) = -^-> llflmll , [/i> 


where ||c|| 2 = YT=i I (a,-, a m ) | 2 . The inequality (71) holds. By 
the definition of / 4 , showing the inequality (72) is equivalent 
to proving 


|c|| 2 < ||a m || 2 ||a m - 1 || 2 . 


(74) 


This is true because 


Proof We first need to show 

II a m II 2 f^Uif/UmW 2 ) 

Clearly, this is true because of (71). Next we need to show 


(76) 


M 

(II &m— 1 

ll 2 /||a m || 2 ] 

1 

max] 

| II &m— 

2 II 2 > / 3 I 

(ik.-illVilaJI 2 )} 


(77) 


+ 1 


O-m -1 I 


I fl-m- ill /II O-m—2 I 
To this end, it is suffice to prove 


p-1 < llflmll 2 + \\a m -i\' 2 
dm 


This is exactly (72). The proof is complete. 


(78) 

n 


The lower bound in (70) can be rewritten in terms of the 
hubgap of A as follows: 


/ 4 (||flil| 2 /||a,„|| 2 ) i 

— -= ;HG r(Q) 


I O-m - 1 I 


1 + 


1 + 


4II ell 


”\ II O-n 


■ (79) 


The upper bound in (70) can be rewritten in terms of the 
hubgap of A as follows: 


M 

^ II Mm- 1 

Il 2 /Ila m ll 2 ) 

1 

max] 

| II &m— 

2 1 2 , / 3 I 

( km-lll 2 / |flm I 2 )} 


/ 4 (||fl m - 1 || 2 /||a m || 2 ) 


I ^m-2 I 


= -(HGi(A) + 1)HG 2 (A) 


llaJI 2 > 


m -1 


s 




2 


> 


m -1 


s 





(75) 


The first inequality of above is from the orthogonality of a ; -, 
j = 1 , ..., m - 1 while the second inequality is from ||fli || < 

11a 2 11 < ■ ■ ■ < ||a m _i ||. This completes the proof. jjjj 


The following result holds. 

Proposition 19. Let Q in ( 66 ) be the arrowhead matrix 
associated with a hub matrix A. Assume that 0 < ||a 4 1| < 
l|a 2 II < ■■■ < Ikm-ill < \\a m \\, where a jt j = l,...,m 
are columns of A. Then the bounds of the EGi(Q) in (70) are 
tighter than those in (69). 


(\\a m \\ 2 -\\a m M\ 2 ) 2 ' 

(80) 

To compare these bounds to Rung and Suter [9], set || c|| 2 = 
0, and the bounds for EigGapj(Q) in (70) become 

HGi(A)<EGi(Q)<HGi(A)HG 2 (A). (81) 


+ ^(HG,(A)-1)HG 2 (A) 

A 


Under these conditions, the lower bound agrees with Rung 
and Suter while the upper bound is tighter. 

Let A E K" xm be an ordered hub matrix. Let A E 
K nxk be a hub matrix obtained by removing the first n - k 
nonhub columns of A with the smallest Euclidean norms. 
This corresponds to the MIMO beamforming problem by 
comparing k of m transmitting antennas with the largest 
signal-to-noise ratio (see [9]). The ratio MQ)/M(Q) with 
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Q = A‘A describes the relative performance of the resulting 
systems. It was shown in [9] that for k > 2 


h( 

112 " WQ) 


I dm II +11 a m -\ II 


(82) 


II Am II +11 1 I 

By Theorems 11 and 12, we have 

f4(\\am-k + i\\ 2 /\\a m \\ 2 ) < h(Q) U(\\cL m -A 2 /\\a m \\ 2 ) 
/ 4 (||fl m -il| 2 /||fl m || 2 ) “'UQ)- / 4 (||fl 1 || 2 /||fl m || 2 ) 

(83) 

By Lemma 18, the lower and upper bounds for the ratio 
Ait(Q)/A m (Q) in (83) are better than those in (82). In 
particular, when k = 1, the matrix Q corresponds to the 
hub, as such, it reduces to Q = [||a m || 2 ]; hence, Ai(Q) = 
||a m || 2 . Therefore, an estimate of the quantity ||fl m || 2 /A m (Q) 
was given in [9] as follows: 

ii 112 ii 112 


l|flml| 2 +l|fl m -ll| 2 A m (Q) 
By Theorems 11 and 12, we have 

II 2 < IlflJl 2 < 


< 1. 


Clyy 


(84) 


(85) 


/ 4 (lk„-ill 2 /lkJ 2 ) / 4 (llcnll 2 /llfl m ll 2 )' 

Again, by Lemma 18, the lower and upper bounds for the 
ratio l|flmll 2 /A m (Q) in (85) are better than those in (84). We 
can simply view (84) and (85) as degenerate forms of (82) 
and (83), respectively. 

5. Numerical Examples 

In this section, we will numerically compare the lower 
and upper bounds of eigenvalues for arrowhead matrices 
estimated by three approaches. The first approach is due to 
Cauchy, denoted by C and based on (2)-(4). The second 
approach, denoted by SS, is based on eigenvalue bounds 
provided by Theorems 11 and 12. The third approach, 
denoted by WS, is based on Wolkowicz-Styan’s lower and 
upper bounds for the largest and smallest eigenvalues of a 
symmetric matrix [15]. These WS bounds are given by 


a - sp < Ai(Q) < a — 
a + ^ < A m (Q) < a + sp, 


( 86 ) 


where Q e 


|mxm 
c 2 


is symmetric, p = \!m- 1, a = 


trace(Q)/m, ands 2 = trace(Q f Q)/m - a 2 . 

Example 20. Consider the directed graph in Figure 1, which 
might be used to represent a MIMO communication scheme. 

The adjacency matrix for the directed graph is 



which is a hub matrix with the right-most column cor¬ 
responding to node 4 as the hub column. The associated 
arrowhead matrix Q = A T A is 


Q = 


10 0 1 
0 2 0 2 
0 0 11 
12 14 


( 88 ) 


The eigenvalues of Q are 0,1,1.4384,5.5616. Corollary 3 
implies 1 being the eigenvalue of Q. By Corollary 4, the 
following matrix 


Q = 


1 0 V2 
0 2 2 
4 2 4 


(89) 


should have eigenvalues of Ai(Q) = 0, Xi{Q) = 1.4384, 
Aj(Q) = 5.5616. The C bounds, SS bounds, and WS bounds 
for the eigenvalues of the matrix Q are listed in Table 1. For 
Ai (Q), the lower SS bound is best; next is the lower C bound 
followed by the lower WS bound; the upper SS bound is best; 
next is the upper WS bound followed by the upper C bound. 
For A 2 (Q), SS bounds and C bounds are the same as the C 
bounds. For A 3 (Q), the lower SS bound is best; the lower C 
bound and WS bound are the same; the upper SS bound is 
best; next is the upper WS bound followed by the upper C 
bound. In conclusion, the SS bounds are best. 

The bounds of the eigengap of Q provided by (69) are 


2 < EGi (Q) < 6 
while the bounds provided by (70) are 

2.6861 < EGi(Q) < 5.6458. 


(90) 


(91) 



0 10 1' 

Example 21. We consider an arrowhea 


0 0 11 


‘ 1 0 V2 

A = 


, (87) 



0 10 1 

Q = 

0 6 2 


10 0 1 


i/2 2 7 


Therefore, the bounds by (70) are tighter than those by (69). 
This numerically justifies Proposition 19. 


( 92 ) 
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Table 1: C bounds, SS bounds, and WS bounds for Example 20. 


Eigenvalues 


C bounds 


SS bounds 



WS bounds 

Lower 


Upper 

Lower 

Upper 

Lower 

Upper 

Ai(Q) = 0 

-0.4142 


1 

-0.3723 

0.3542 

-1 

0.6667 

A 2 (Q) = 1.4384 

1 


2 

1 

2 



A 3 (Q) = 5.5616 

4 


7.4142 

5.3723 

5.6458 

4 

5.6667 

Table 2: C bounds, SS bounds, and WS bounds for Example 21. 

Eigenvalues 


C bounds 


SS bounds 



WS bounds 

Lower 


Upper 

Lower 

Upper 

Lower 

Upper 

A i(Q) = 0.6435 

-0.4142 


1 

0.1270 

1 

0 

2.3333 

A 2 (Q) = 4.6304 

1 


6 

4 

6 



A 3 (Q) = 8.7261 

6 


10.4142 

7.8730 

9 

7 

9.3333 


Table 3: C bounds, SS bounds, and WS bounds for Example 22. 


Eigenvalues 

Lower 

C bounds 

Upper 

Lower 

SS bounds 

Upper 

WS bounds 

Lower 

Upper 

Ai(Q) = 0.6192 

0.5000 

1 

0.2753 

0.8820 

0.2137 

0.9402 

A 2 (Q) = 1.3183 

1 

2 

1 

2 



A 3 (Q) = 3.0624 

2 

3.5000 

2.7247 

3.1180 

2.3931 

3.1196 


The eigenvalues of Q and the corresponding C bounds, 
SS bounds, and WS bounds for the eigenvalues of Q are listed 
in Table 2. For Ai (Q), the lower SS bound is the best; next is 
the lower WS bound followed by the lower C bound; the SS 
upper bound and the C upper bound are the same and are 
better than the upper WS bound. For the lower SS 

bound is better than the lower C bound; the upper SS bound 
is the same as the upper C bound. For A 3 (Q), the lower SS 
bound is the best; next is the lower WS bound followed by 
the lower C bound; the upper SS bound is the best; next is 
the upper WS bound followed by the upper C bound. 

Example 22. We consider an arrowhead matrix Q as follows: 


Q 



0 2 1 


(93) 


6. Conclusions 

Motivated by the need to more accurately estimate eigengaps 
of the system matrices associated with hub matrices, this 
paper provides an efficient way to estimate the lower and 
upper bounds of arrowhead matrices. Improved lower and 
upper bounds for the eigengaps of the system matrices are 
developed. We applied these results to a wireless computation 
application, and subsequently we presented several numeri¬ 
cal examples. In the future, we will plan to extend our results 
to hub dominant matrices, which will allow hub matrices 
with correlated nonhub columns. 
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